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Abstract. It has been suggested that underwater 
neutrino telescopes could detect muons from gamma 
ray showers. Antares' ability to detect high energy 
muons produced by TeV photons is discussed in the 
light of a full Monte Carlo study. It is shown that 
currently known sources would be hardly detectable. 
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I. Introduction 

The last decade has been fruitful in terms of high 
energy astronomy. More than 80 gamma rays sources 
have been found to emit in the TeV range [1], thanks 
to Imaging Atmospheric Cerenkov Telescopes (IACTs) 
such as Cangaroo, HEGRA, H.E.S.S, Magic or 
Veritas /Whipple 0. 

There have also been many discussions about the 
possibility to detect muons produced by high energy 
gamma rays in underground, underice or underwater 
neutrino telescopes |3|[4||6|. In contrast to upward- 
going muons from neutrinos, which are the primary 
purpose of such a telescope, downward-going muons 
from gamma rays suffer from a high atmospheric muon 
background. Therefore the sensitivity of a neutrino tele- 
scope to gamma ray induced muons is quite lower than 
IACTs'. However, it has the advantage of monitoring 
continuously all directions. In addition to their physics 
potential, muons from gamma rays may also offer 
calibration benefits in terms of pointing accuracy and 
angular resolution. 

Gamma ray showers are believed to be muon poor, 
but there are at least three processes by which a photon 
can produce muons : photoproduction, muon pair pro- 
duction and charm decay. The first process involves the 
(semi)leptonic decay of a pion or a kaon produced by the 
interaction of the photon with an atmospheric nucleus. 
Such muons are said to be conventional. The second 
process is self-explanatory, and its final particles are 
referred to as direct muons. The final case corresponds 
to the (semi)leptonic decay of a photoproduced charm 
meson, and secondary muons are called prompt muons. 
The prompt muon production was not taken into account 
in this work since it was not implemented in the soft- 
ware used for the Monte Carlo production. The charm 
production involves QCD processes that are not fully 
understood, but measurements at HERA have shown that 
at photon energies of several TeV charm production is 
significant |5|. 



Some calculations have estimated that the muon flux 
from gamma ray sources could be sufficient for neutrino 
telescopes to detect them. However, most attempts to 
estimate this muon flux rely on one-dimensional analytic 
models, and do not take into account the muon prop- 
agation from sea level to the detector and the detector 
sensitivity. A first attempt to estimate the underwater flux 
using a Monte Carlo simulation, without considering 
detection efficiency, has found gamma ray sources to be 
hardly detectable by a neutrino telescope (6). 

In this paper, a full Monte Carlo simulation, including 
Cerenkov light detection in realistic background condi- 
tions and track reconstruction, is presented, within the 
Antares framework. The expected number of events 
from the main sources of interest are presented. 

II. The Antares detector 

The Mediterranean sea currently houses the first oper- 
ational undersea neutrino telescope, and also the largest 
neutrino telescope in the Northern hemisphere, namely 
Antares (Astronomy with a Neutrino Telescope and 
Abyss environmental RESearch) Q. Its full configura- 
tion has been completed in May 2008, though data has 
been taken with partial detector configurations since the 
first line was in water, in March 2006. 

Antares' main focus is to detect astrophysical neu- 
trinos, thanks to the Cerenkov light produced in water by 
muons resulting from the interactions of neutrinos with 
the Earth. Because of the atmospheric muon background, 
Antares field of view is the Southern hemisphere, and 
in particular the Galactic center. 

Installed at 40 km off Toulon, in France (40°50'N, 
6° 10' E), Antares comprises 12 vertical detection lines 
positioned on the Mediterranean sea bed, at about 
2500 m depth. Each line hosts up to 25 floors of three 
10" photomultiplier tubes (PMTs) regularly distributed 
over 350 m, the lowest floor being 100 m above the 
sea bed. On a given floor, the three PMTs are placed 
symmetrically around the vertical string and look down- 
wards at 45° in order to optimize the collection of 
Cerenkov photons from upgoing muons rather than from 
downgoing muons lf8~l . 

The lines are separated from each other by approxi- 
mately 70 m, and set as a regular octagon surrounding 
a square. An instrumented line intended to monitor the 
environmental conditions completes the apparatus. 

The sea current induced displacements of the lines 
with regard to their vertical axis do not exceed a few 
meters, and are monitored in real time using compasses, 
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tiltmeters and hydrophones hosted on each line. A posi- 
tion accuracy of about 10 cm for each PMT is obtained. 

An electro-optical cable transfers the electronic read- 
out of the whole detector to shore, where digitized 
informations are processed in a computer farm. 

III. Monte Carlo simulation 

Extensive Air Shower have been simulated using 
Corsika v6.720 [9|. High energy hadronic interactions 
are modeled through QGSJET01 [ 10 1, while electro- 
magnetic interactions are processed through EGS4 ifTTI . 
QGSJET01 is found to be the most conservative model 
in comparison to SIBYLL, VENUS and QGSJET-II ifTOl . 
regarding the number of photons creating high energy 
muons (using VENUS leads to a 7% rise, assuming a 
E" 1 flux, in the [l;100]TeV range). However, the effect 
of this increase at the depth of ANTARES still has to be 
investigated, the muon range being energy dependent. 
The energy range considered in the present work goes 
from ITeV to lOOTeV. 

MUSIC has been used for the propagation of muons 
in water 0/2). The ANTARES Monte Carlo simulation 
chain then allows for simulation of Cerenkov light in 
the detector, taking into account, in particular, the water 
properties and the PMTs angular acceptance lfl3l . It 
also allows for the addition of realistic bioluminescence 
background using real data streams. 

In this work the data used for the bioluminescence 
background corresponds to golden running conditions : 
runs are selected where the baseline of raw counting 
rates and the fraction of burst|3 are low (about 60 kHz 
and less than 20%, respectively). 

Finally, the events are reconstructed using ANTARES 
standard reconstruction strategy. It has to be noticed that 
this strategy is optimized for upgoing events. The results 
presented here might thus be slightly enhanced using a 
dedicated strategy. On the other hand, the cut made on 
the reconstruction quality is very loose, so the effect of 
hardening the quality cut may compensate the effect of 
improving the reconstruction strategy. 

IV. Sources of interest 

In order to have a reasonable probability to reach the 
depth of the ANTARES detector, a downgoing vertical 
muon must be more energetic than 1 TeV at sea level : 
the muon probability to survive to a 2200 m depth in wa- 
ter is less than 70% for a 1 TeV muon (13% at 700 GeV). 
Hence only TeV gamma ray sources may be seen by 
ANTARES. More than 80 gamma ray sources have been 
detected in the TeV range by IACTs JT]. However, not 
all of them are good candidates for ANTARES : most of 
them are located in the galactic plane, which is not in 
ANTARES field of view|3. Moreover, weak and/or soft 
fluxes are not likely to produce enough muons. 

'The burst fraction is defined as the ratio between the time when 
the counting rate is higher than 250 kHz and the overall time. 

2 The gamma rays field of view being the Northern hemisphere, as 
opposed to the neutrinos field of view. 



Fortunately, several of the most powerful sources are 
visible by ANTARES, including the so-called "standard 
candle", the Crab pulsar. Characteristics of the most 
interesting candidates in terms of fluxes and visibility 
are summarized in table U In addition to the Crab, three 
extragalactic sources have been selected. 

Though most of these sources are variable or flaring 
sources, they are known to have long periods of high 
activity, which make them more promising over a long 
period than most steady sources lfT4ll lfT5ll 161. 



source 


visibility 


mean zenith 


type 


Crab 


62% 


51.7 


PWN 


1ES 1959+650 


100% 


49.7 


HBL 


Mkn 501 


78% 


49.4 


HBL 


Mkn 421 


76% 


49.2 


HBL 



TABLE I 

Characteristics of ANTARES best gamma ray sources. HBL stands for 
High frequency peaked BL Lac object, and PWN for Pulsar Wind 
Nebula. 



V. Simulation results 

The number of detected and reconstructed events 
depends on several parameters, such as the precise level 
of background, the trigger strategy, the reconstruction 
strategy and the source flux parametrization. Therefore 
only range estimates are given. They are reported in 
table HU assuming a 100% visibility over one year. Most 
excentric parametrizations have been omitted. 

It is found that only a few photons can be expected to 
be seen by ANTARES : less than ten events per year are 
reconstructed for the Crab in realistic conditions, though 
a few tens produce hits on the detector. 

In comparison, a rough estimate on data with sim- 
ilar bioluminescence conditions gives 1.1 x 10 5 (resp. 
4.1 x 10 4 ) reconstructed background events (atmospheric 
muons) within a one degree cone of the 10 degrees 
zenith angle direction (resp. 40 degrees), and 2.1 x 10 6 
(resp. 7.3 x 10 5 ) background events within a 5 degrees 
cone. 

It seems thus not reasonable to expect ANTARES to 
extract any gamma signal from the background under 
these conditions for any known flux. Though Markarian 
501 may seem promising, the upper limit actually corre- 
sponds to high state fluxes parametrizations on dayscale 
variations [16). A more precise selection of the fluxes 
and a study of the significance of the expected number 
of muons are still to be done. 

However, these estimates are not so bad as they 
seem. First, the simulation is conservative in terms of 
photoproduction cross-section [17| and does not take 
muon production from charm decay into account. In 
addition, background discrimination has not yet been 
investigated. In particular, the muon poorness of gamma 
ray showers may help to reduce the atmospheric back- 
ground : by rejecting multimuon events, one can improve 
both the signal to noise ratio and the angular resolution. 
If achievable, the muon pair tagging may also improve 
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source 


patm 


F= ca (xl0- 3 ) 


N det 


N dot 
iu-t- 


^rcco 


Crab 


5-8 


0.2-0.8 


30-70 


20-45 


1-4 






0.4-1 


20-50 


15-35 


1-3 


1ES 1959+650 


0.8-30 


0.1-2 


3-150 


2-100 


0.2-8 






0.1-2 


2-100 


1-70 


0.1-5 


IV ^ /IT 1 

MKn 421 


1.5-45 


0. 1-4 


5-330 


3-230 


0.2-20 






0.1-5 


2-260 


2-190 


0.1-15 


Mkn 501 


1.5-40 


0.1-15 


6-1350 


4-950 


0.3-90 






0.1-15 


3-1200 


2-880 


0.1-80 



TABLE II 

Number of photons/muons produced by several gamma ray sources at different levels, during one year, assuming a 100% visibility, for primaries 
in the [1;100] TeV energy range. When relevant, straight font corresponds to a 10 degrees zenith angle, while italic stands for a 40 degrees zenith 
angle. Fluxes of photons at the top of the atmosphere (F atm j and at sea level (F° ca j are expressed in m~ 2 . N^* is the number of photons which 
produce more than X hits on the detector PMTs, and N rcco corresponds to the number of reconstructed events in realistic bioluminescence 

conditions. 



the background rejection. Moreover, a dedicated recon- 
struction strategy could increase the number of detected 
photons. Finally, galactic sources such as the Crab are 
not subject to the universe opacity above 100 TeV, and 
the extension of their spectra to higher energy could 
lead to reasonable numbers of detectable photons. Short 
and powerful bursts are not to be excluded either, 
the associated background being in such cases almost 
negligible. 

VI. Conclusion 

A complete Monte Carlo study has been processed in 
order to estimate Antares' ability to detect downgoing 
muons from gamma ray sources. 

It has been found that ANTARES is not likely to detect 
any of the currently known sources, unless they show 
some unexpected behaviour. However the conservative 
estimates computed in this work show that the gamma 
ray astronomy field is not completely out of reach of 
underwater neutrino telescopes, at least for the next 
generation of detectors. 

This study will be refined and extended to the km 3 - 
scale successor of ANTARES, namely KM3NeT, which 
is currently being designed ifTSl . 
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